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Abstract. It has been shown in several recent papers that the six Doppler 
data streams obtained from a triangular LISA configuration can be combined by 
appropriately delaying the data streams for cancelling the laser frequency noise. Raw 
laser noise is several orders of magnitude above the other noises and thus it is essential 
to bring it down to the level of other noises such as shot, acceleration, etc. A rigorous 
and systematic formalism using the powerful techniques of computational commutative 
algebra was developed which generates in principle all the data combinations cancelling 
the laser frequency noise. The relevant data combinations form a first module of 
syzygies. 

In this paper we use this formalism to advantage for optimising the sensitivity of 
LISA by analysing the noise and signal covariance matrices. The signal covariance 
matrix is calculated for binaries whose frequency changes at most adiabatically and 
the signal is averaged over polarisations and directions. We then present the extremal 
SNR curves for all the data combinations in the module. They correspond to the 
eigenvectors of the noise and signal covariance matrices. A LISA 'network' SNR is also 
computed by combining the outputs of the eigenvectors. We show that substantial 
gains in sensitivity can be obtained by employing these strategies. The maximum 
SNR curve can yield an improvement upto 70 % over the Michclson, mainly at high 
frequencies, while the improvement using the network SNR ranges from 40 % to over 
100 %. 

Finally, we describe a simple toy model, in which LISA rotates in a plane. In this 
analysis, we estimate the improvement in the LISA sensitivity, if one switches from 
one data combination to another as it rotates. Here the improvement in sensitivity, if 
one switches optimally over three cyclic data combinations of the eigenvector is about 
55 % on an average over the LISA band-width. The corresponding SNR improvement 
increases to 60 %, if one maximises over the module. 
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1. Introduction 

The proposed space-based mission of gravitational wave (GW) detection - the Laser 
Interferometric Space Antenna (LISA) - consists of three identical spacecrafts forming an 
equilateral triangle of side 5 x 10 6 km following heliocentric orbits. LISA is thus a giant 
interferometric configuration with three arms which will give independent information 
on GW polarisations. The implementation of this huge Michelson interferometer is 
quite different from the ground based interferometers like LIGO, VIRGO etc. Unlike in 
the ground based interferometric detectors, it is not feasible to bounce the laser beams 
because of the large armlengths involved. A special Doppler tracking scheme is used 
to track the space crafts with laser beams. This exchange of laser beams between the 
three space-craft result in six Doppler data streams. 

LISA sensitivity is limited by several noise sources. A major noise source is the 
laser phase noise which arises due to phase fluctuations of the master laser. Amongst the 
important noise sources, laser phase noise is expected to be several orders of magnitude 
larger than other noises in the instrument. The current stabilisation schemes estimate 
(this estimate may improve in the future) this noise to about Au/vq ~ 1CT 13 /\/Hz, 
where z/ is the frequency of the laser and Au the fluctuation in frequency. If the laser 
frequency noise can be suppressed then the noise floor is determined by the optical- 
path noise which fakes fluctuations in the lengths of optical paths and the residual 
acceleration of proof masses resulting from imperfect shielding of the drag-free system. 
The noise floor is then at an effective GW strain sensitivity h ~ 1CT 21 or 1CT 22 . Thus, 
cancelling the laser frequency noise is vital if LISA is to reach the requisite sensitivity. 
Since it is impossible to maintain equal distances between space-craft, cancellation of 
laser frequency noise is a non-trivial problem. Several schemes have been proposed 
to combat this noise. In these schemes |l|, |J, the data streams are combined with 
appropriate time delays in order to cancel the laser frequency noise. In our earlier 
work ||, henceforth referred to as paper I, we had presented a systematic and rigorous 
method using commutative algebra which generates all the data combinations cancelling 
the laser frequency noise. These data combinations form a module over a polynomial 
ring, well known in the literature, as the first module of syzygies. We obtained the 
generators of this module and hence the entire set of relevant data combinations could 
be generated conveniently. More importantly, we note that this method is general and 
can be extended to space-missions with more than three spacecrafts. 

In this paper we employ our previously set up formalism for two important 
applications: We compute the noise covariance matrix for LISA and obtain its 
eigenvectors and eigenvalues. The signal covariance matrix is computed for binaries 
whose frequency changes at most adiabatically (the monochromatic case is included) 
and for which the signal is averaged over polarisations and directions. Here adiabatic 
means that the signal response, the noise and hence the SNR change imperceptibly 
even if the GW source changes frequency during the observation time. Thus, even 
though the results are presented at each fixed frequency, the sources need not be strictly 
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monochromatic, and apply to a wider class of sources. The signal covariance matrix has 
the same eigenvectors as the noise covariance matrix which results in computational 
simplification. We show that the signal-to-noise (SNR) for any data combination in 
the module, then lies between an upper and a lower bound. The upper and lower 
bounds of the SNR are functions of frequency which are just the SNR curves of the 
eigenvectors. The extremisation - both maximisation and minimisation - of SNR is 
important for different purposes; maximisation is important for the detection and 
parameter estimation of a GW source, while minimisation is important for the purpose 
of distinguishing the GW confusion noise from the instrumental noise Q. We further 
show that the bounding SNR curves have multiple intersections within the band-width 
of LISA (10~ 4 — 1) Hz. The improvement of SNR of the upper-bound over the Michelson 
combination goes upto 70 % , but only at high frequencies > 5 mHz. At low frequencies 
< 5 mHz, both have the same sensitivity. Since the eigenvectors are independent 
random variables, a 'network' SNR of independent detectors || can be constructed 
from the likelihood considerations which gives an improvement between a/2 and \/3 
over the maximum of SNRs of the eigenvectors. The improvement over the Michelson 
combination is about 40 % at low frequencies < 3 mHz and rises above 100 % at high 
frequencies. We may note that, some of our results are in agreement with independent 
and simultaneous calculations by Prince et al fl6|. 

Tracking a GW source fixed on the celestial sphere for determining the SNR 
constitutes a non-trivial problem, which will require substantial analysis when writing 
codes. Here, we consider a simple toy model of LISA rotating in its own plane. Our 
goal is to obtain a rough estimate in the improvement of SNR by optimally switching 
from one data combination to another. We show that, if we implement this strategy, it 
is possible to improve the SNR by about 55 %. 

2. The module of syzygies of time delayed data combinations 

In this section, we briefly summarise the main results of paper I. It has been shown 
earlier in the literature § || , that the laser phase noise and the optical bench motion 
noise in the LISA can be suppressed by combining the six beams appropriately delayed 
across arms of the interferometer. The six data streams are labelled as U l and V" 1 , 
i = 1,2, 3. The geometry of the LISA along with six laser beams is shown in the Figure 
-|l]. The beams U l go clockwise while the V 1 go counter-clockwise. There are also six 
other beams connecting the two optical benches on each space-craft (which are not 
shown in the figure). In the SNR analysis only the difference between these beams is 
important. These differences are denoted by Z l where the index i corresponds to each 
space-craft labelled i. 

In paper I, we have shown that all the data combinations which cancel the laser 
frequency noise and the bench motion noise form a module over a ring of polynomials 
in the three time delay operators Ei, where E iy i = 1,2,3 represent delay operators of 
the light travel time along the three arms with lengths Lj. Thus for the data stream 
x(t) : Eix{t) = x{t — Li) (the speed of light is set to unity); the operator Ei delays the 
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Figure 1. The LISA constellation consist of three spacecraft each carrying two optical 
bench systems i and i* , i — 1, 2, 3. They exchange six laser beams which are represented 
by Wand V 1 in the figure. The noise cancellation data combinations are obtained by 
appropriately delaying these six beams across the length of arms Li. 




data stream by amount Lj. For the LISA configuration Lj ~ 16.7 seconds, corresponding 
to an arm- length of 5 million km. An important advantage of this formalism is that, 
one can generate the entire module from the generators; linear combinations of the 
generators with polynomial coefficients in the ring generate this noise cancellation 
module. Several sets of generators have been listed in paper I. Based on the physical 
considerations or for the purpose of computational convenience, we may choose one set 
over another. For the purpose of extremising the SNR, we choose the set of generators: 
a, /3,7 and ( (notation followed from |l], 0, [3|]). The a, /3, 7 are cyclic permutations of 
each other. This symmetry comes in useful when computing scalar products between 
them and also in diagonalising the noise, and signal covariance matrices defined in the 
next section. Following paper I, we list this generating set as 9-tuples of polynomials 
(pi, q iy 7™j). The polynomials (pi, qi, r^) are polynomials in the variables Ei and act on the 
data streams (V*, U l , Z l ) respectively. The linear combination piV 1 + qiU l +ViZ l cancels 
the laser frequency noise and optical bench motion noise when the 9-tuple (pi, q^, rj) is 
an element of the module of syzygies. The generators of the module are given by: 

a = (1, Es,EiEs, 1, E1E2, E2, —(1 + E1E2E3), —(E1E2 + E3), —(E1E3 + £2)) , 
P = (E1E2, 1, Ex,Es, 1, E2E3, —(E1E2 + E3), —(1 + E1E2E3), —(Ei + E2E3)) , 
7 = (E 2 ,E 2 E 3 , 1, E 3 E 1: Ex, 1, -{E 2 + £i# 3 ), -(#1 + E 2 E 3 ), -(1 + E 1 E 2 E 3 )) , 
C = {Ei,E2, E% , Ei , E2 , E3 , —(Ei + E2E3), —(E 2 + E1E3), —(E3 + E1E2)) . (1) 

We exhibit another set of generators { } which was obtained from the software 
package Macaulay 2 [^]. They can be related to the above set of generators by: 
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x^ = C, 

X^ = (3. (2) 

Note that these sets of generators are not linearly independent. In particular, the 
set of generators {a,f3, 7, C} [0] obey the following condition : 

(1 - £i£ 2 £ 3 )C = {Ei - E 2 E 3 )a + (E 2 - E 1 E 3 )/3 + (E 3 - E X E 2 ) 1 . (3) 
When maximising the SNR in the Fourier space, this relation allows us to eliminate one 
of the generators at almost all frequencies, except when the product E\E 2 E 3 = 1. Note 
that ExE 2 E 3 is just the total time-delay L\ + L 2 + L 3 around the LISA triangle. For 
the purposes of this paper, we assume that all the arms of LISA are of equal length i. e 
L\ = L 2 = L 3 = L. In the Fourier domain, i.e, E\ — E 2 — E 3 — E — e lQL and the 
operator polynomials become actual polynomials. One can then solve for ( in terms 
of a,/3, 7, except at the frequencies Q, when e 3lUL = 1. Taking L ~ 5 x 10 5 km, the 
smallest such frequency / = Q/2tt is ~ 20 mHz. Thus, 

^TT§Te> + + -^ (4) 

and can be effectively eliminated while extremising SNR, except at the singular 
frequencies. Since SNR computation can be successfully carried out for frequencies 
arbitrarily close to the singular frequencies, for the computation, the singularities do 
not seem to be important. In the analysis that follows, we use only the three generators 
{a, /3, 7}- 



3. Strategies for improving the effective sensitivity of LISA 

In this section, we show that the set of generators {a, j3, 7} can be combined into a 
new set, consisting of 'orthogonal' eigenvectors. The noise covariance matrix naturally 
defines a positive definite, non-degenerate bilinear form, which serves as a scalar 
product or a metric. Orthogonality between eigenvectors is defined in terms of this 
metric. Physically this means that the noises of the eigenvectors are uncorrelated 
with each other. The eigenvectors are easily obtained by diagonalising the noise 
covariance matrix. The averaged signal matrix that we consider here has the same 
form as the noise covariance matrix and consequently has the same eigenvectors. Thus 
this set of eigenvectors simultaneously diagonalises both matrices constituting the 
SNR and simplifies the analysis that follows. An important observation here is that 
the eigenvectors are independent observables. They represent therefore statistically 
independent detectors (so far as instrumental noises are concerned), and they can 
be treated as a network of detectors. Furthermore, they can be combined in a root 
mean square fashion to yield a 'detector network statistic' H to yield a much improved 
sensitivity. 
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3.1. The noise covariance matrix 



Following the formalism in paper I, we define noise vectors in the Fourier domain 
N"\ I = 1, 2, 3 for each of the generators {a, (3, 7} respectively, over the 12 dimensional 
complex space C 12 , 



S sh p { -\^S sh q, 



(J) 



(5) 



where S p f(f) and S opt (f) are power spectral densities (psd) of the proof mass residual 
motion and the optical path noise respectively. The polynomials (p^ , gj ^ , r^ ) , (now 
actual polynomials in the Fourier domain) corresponding to the generators a, j3 and 
7 are given in the equation (0). We take S pf (f) = 2.5 x KT 48 [f/1 Hz]' 2 Hz" 1 and 
S°p\f) = 1.8 x 1(T 37 [f/1 Hz} 2 Hz- 1 following the literature §. It can be easily shown 
that for a given data combination, the norm of the noise vector represents its noise psd. 
The noise covariance matrix for the generators {a, /3, 7} is defined as J\f( IJ ^ = N^-N^ J '* 
and takes the simple form, 



rid 

Tin 



n 
n 
rid 



(6) 



We note that because of the cyclic symmetry, the diagonal elements ■ N^* are 
equal to each other - denoted by rid- Similarly, all the off-diagonal elements ■ iV^*, 
for I J are also equal to each other and which we denote by n D . This was the 
reason a generating set possessing symmetry properties was chosen in the first place. 
A matrix with this form has two degenerate eigenvalues. Thus, the eigenvalues of the 
noise covariance matrix are given by, 

m = n 2 = n d — n Q and n 3 = n d + 2n a . (7) 

Since two of the eigenvalues are degenerate we need to systematically adopt a procedure 
for choosing the linearly independent and orthonormal set of eigenvectors. This choice 
is not unique. One such choice gives the following matrix M. 

11 [2 



M 



Vg Ve V 3 
__L 1 n 

111 
V3 V3 \/3 



which diagonalises the noise matrix J\f( IJ \ that is Ai 
eigenvalues as diagonal elements. The eigenvectors are: 

' 0-27), 



(8) 



J\f ■ M. 1 is diagonal, with 



y(i) 
y(2) 



1 

V2 



[a 



a) 



= — ( a + (3 + 1 ) . 



(9) 



We find that the data combination is proportional to the symmetric Sagnac 
combination ( and has the same SNR as that of £. 
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3.2. The signal covariance matrix 

The response of a GW signal for a given data combination is computed in paper I. The 
response is conveniently expressed in the Fourier domain and is given by, 

h(Q) = £ [p t (F Vi . + h + + F v . x h x ) + q t (F v . + h + + Fu. x h x )} (fi) . (10) 

i=l 

Here, Fy.. +/x and Fu i . +/X are the antenna pattern functions. We note that the signal 



»;+/x 



depends only on the first six entries (p i: of the 9-tuple describing a data combination. 
So while dealing with the signal response, we only consider the 6-tuple P = (pi,qi) of 
the 9-tuple describing a data combination. We apply this formalism to a binary source 
which may be adiabatically changing in frequency. The two GW amplitudes of such a 
source at frequency Q are given by, 

, , (l + cos 2 e n t . . n ,\ 

n + {U) = A I cos 2ip — i cos e siri2%p I , 

(1 + cos 2 e \ 
sin2tp — i cos e cos 2ip j . 

Here, the parameters e and ip describe the orientation of the source and enter into the 
expressions for the polarisation amplitudes. The direction of the source on the celestial 
sphere is given by the angles 9 and 0. In order to organise the calculations, we also 
define the detector response 6-tuple as, 

R = (F Vi;+ h + + F v . x h x , F Ut;+ h + + Fu. x h x ) . (11) 
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Both P and R will be considered as row vectors for the purposes of defining matrix 
products. In order to analyse the signal covariance matrix, it is useful to define a scalar 
product. For two data combinations P and Q (considered as row vectors), we define the 
scalar product as follows: 

<P,Q> = P-IZ-Q], (12) 

where, TZ = PJ ■ R is a Hermitian matrix of detector responses and the 'dot' denotes 
the matrix product. The norm of the vector P is then given by, 

\\P\\ 2 = <P,P> . (13) 

The norm of P is the GW response for the data combination described by P. 

The signal covariance matrix for any generating set X' 1 ' (and corresponding P^) 
is then defined as, 

H (U) = {h^ J >) e ^ = (« P (I) , P {J) ») #(¥ , (14) 

where, = PW • PJ and the bracket ( ) e ^e(f> represents the average over the 

polarisations and directions. Taking a, f3 and 7 as the generators, the cyclic symmetry 
between them gives rise to a signal covariance matrix H,( IJ ^ which has the same form as 
the noise covariance matrix J\f( IJS) given in equation (|(J). In this case the nd and n Q are 
replaced by = (h^h^*) and h Q = (h^h^*) respectively. Thus Ti. is diagonalised by 
the similarity transformation A4 and has the same eigenvectors Y^\ The eigenvalues 
of are given by 

hi — h<i — hd — h and h 3 = hj + 2h . (15) 

This simultaneous diagonalisation of both signal and noise matrices is important 
from the point of extremisation of SNR. This forms the content of the next subsection. 
However, we may note that in the formalism developed by Prince et al |J, the 
optimisation is performed without averaging over the source directions and polarisations, 
which results in the GW source matrix of rank 1. Since, the source directions are not 
known in general, we average over these parameters which results in a signal matrix of 
rank 3. 



3. 3. Extremisation of SNR 

An arbitrary data combination can be written as Y = amY^ 1 ', where am are 
polynomials in E. The SNR for this combination is given by, 

Siw= M7TW (16) 

where, ki = \ai\ 2 + |ct 2 1 2 and k 3 = \as\ 2 , and the ranges of ki, fc 3 are from to 00. The 
simplicity of this expression is because of using the new set of the orthogonal generators 
Y^\ It enables us to easily solve the extremisation problem. The maximum and the 
minimum of this expression is governed by the quantity rj where, 

VU) SNRf 3) h 3 /n 3 > [U) 



Improving the Sensitivity of LISA 



9 



Figure 3. Log Log plot of sensitivity S, curve as function of / after averaging over 
polarisation and source directions for a observation period of one year and SNR =5. 
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where SNR (I) denotes the SNR of When rj(f) ^ 1, the SNR is a strictly monotonic 
function of the ratio ki/k 3 . If r](f) > 1 then the SNR of the generator or Y^ 
is greater than that of Y^ 3 \ Then for these frequencies Y^ or Y^ (or any linear 
combination of them) yields the maximum SNR while y (3) yields minimum SNR. When 
r](f) < 1, the opposite is true: Y^ or Y^ (or any linear combination of them) yield the 
minimum SNR while Y^ yields maximum SNR. The remaining case, when rj(f) = 1 
all the 7^ have the same SNR. Since the extremum values of SNR are only attained 
by the eigenvectors, the corresponding SNR curves constitute the bounding curves for 
any linear combination of Y^ s. So our results determine the limiting sensitivities of 
data combinations cancelling laser frequency noise and optical bench motion noise. It is 
interesting to note that at the frequencies where the bounding curves intersect, all the 
data combinations belonging to the module have the same SNR. 

In the lower frequency range (/ < l5mHz) 7 the combination F^(same as that 
of the combination () has a very low sensitivity to the gravitational wave signal and 
the generators Y^ (or Y^) has maximum sensitivity to the signal. However, at 
high frequencies, the sensitivity curves of Y^ and Y^ intersect at several frequencies 
eg. ~ 27 mHz, 39 mHz etc. While computing l^S 1 ^ and the average sensitivity, we 
assume a uniform distribution of sources over polarisations and source directions in the 
sky. Averaging over the polarisations is performed analytically and the averaging over 
source directions performed using the Monte Carlo method. The sensitivity, S is defined 
following the reference |§ , S — 5 J g§^, here, B = ~ and T is the observation time which 
we take to be one year. In the Figure |^, we show the plots of S for the basis elements 
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Figure 4. Plot of parameter r\ as a function of /. The points r\ — 1 correspond to the 
frequencies at which all the data combinations have the same SNR. rj >1 corresponds 
to the region in which data combination Y^and Y^> are more sensitive than 
and vice versa. 
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yW, for comparison, we also plot the sensitivity for the Michelson combination X. In 
the Figure |], we plot the ratio rj as a function of the frequency /. The points at which 
rj intersects the line y = 1 corresponds to the points where all the data combinations in 
the module have the same value of SNR. 



3.4.. Operating LISA in a network mode 

In the preceding sections, we have shown that either Y*W, or maximise the SNR 
and they are orthogonal i.e. they are independent random variables. The sensitivity 
of LISA can be further improved because, each of these generators can be realized as 
independent gravitational wave detectors. Here we obtain the network SNR by taking 
Y^ as independent outputs of a network of three detectors. We assume that the 
underlying noise is Gaussian || and the Y^ 1 ' follow the standard normal distribution. 

As shown in ||, if the noise of the individual detectors is uncorrelated then the 
network likelihood ratio is just the product of individual likelihood ratios; the log 
network likelihood ratio is the sum of the individual log likelihood ratios. Moreover, if 
the noise in the individual detectors is Gaussian, the surrogate statistic of the network, 
yields the network SNR as, 



SNR 2 



network 



7=1 



;is) 
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Figure 5. Plots showing the relative improvements (ratios) of SNRs for the three 
cases: (i) Network SNR over the Michelson data combination (solid line), (ii) Network 
SNR over the maximum of Max [Y^'Y^] (dotted line), (iii) Max [yW>y( 3 )] over 
the Michelson (dashed line). Here Max [yW>y( 3 )] is the maximum of the SNR of 
and y( 3 ) over the bandwidth of LISA. 



Network/Michelson 

2.4 - Max[Y ,Y ]/Michelsqri 

Network/Max[Y\Y J ] 




1 \J Wv W 



0.8 

0.0001 0.001 0.01 0.1 

Frequency 



The corresponding sensitivities are shown in the Figure At low frequencies / < 15 
mHz, the improvement of the network SNR over the maximum of Y^ 1 ' is slightly greater 
than v^2- This is because at low frequencies the data combination is not very 
sensitive in comparison with The best improvement in the relative SNR is achieved 
at frequencies where all the data combinations are equally sensitive, that is, when r\ = 1. 
A gain of a factor of \/3 is achieved at these points. In the Figure |^, we have plotted the 
relative improvements in the network SNR with respect to the Michelson combination 
and the maximum of and F (3) . 

4. SNR maximisation for a toy model of the LISA triangle rotating in a 
plane 

The LISA configuration consists of three spacecraft which orbit with complex 
trajectories in order to maintain a stable triangular configuration in an heliocentric 
orbit ||. We note that the transfer functions of the beams U l , V 1 are computed in the 
frame of LISA. Hence, a given GW source which is fixed on the celestial sphere will 
appear to follow a complex trajectory in this frame. This motion is the superposition of 
two rotations. One is just a rotation of the LISA triangle in its own plane. The other is 
its motion around the Sun in which the configuration follows the earth in an earthlike 
orbit. The period of the two rotations is identical and is one year in duration. Since 
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LISA has a non-uniform directional response, it is a non-trivial problem to track the 
apparent motion of a 'fixed' GW source in the LISA frame and then compute its SNR. 

We consider a toy model in which only one motion is taken into account, which is 
the rotation of LISA about an axis orthogonal to the plane containing the three space- 
craft. We ignore the other motion of LISA around the Sun for simplicity. Following 
paper I, we choose the X — Y plane to coincide with the plane formed by the three 
space-craft constituting the LISA triangle. The origin is chosen at the centre of the 
equilateral triangle (for computing the response, we take all the arms to be of equal 
length). The positive X-axis passes through space-craft 2. Figure 1 shows how the 
X — Y axes have been chosen. The Z-axis is orthogonal to this plane with the positive 
direction given by the right-handed convention. Under these assumptions, the source 
would appear to move in the sky in a circular trajectory about the Z-axis with a period 
of one year. In terms of polar coordinates (9, 0), the motion of such a source is uniform 
along with a constant value of 9. If one tracks a source with a single data combination, 
say Y^ l \ because the directional sensitivity of Y^ is non-uniform, it does not track the 
source optimally. Figure (j^j displays the 3-dimensional sensitivity plots of Y^ and Y < - 2 - ) 
as a function of the angles (9, 0) for a monochromatic source at the GW frequency of 1 
niHz and the signal uniformly averaged over the polarisations (this average is computed 
analytically). It is obvious from these plots that the sensitivity is a highly non-uniform 
function of source direction. We consider here the case when the source with GW 
frequency 1 mHz lies in the plane of the LISA triangle, that is, the GW source lies in 
the 9 = | plane. We choose this value because, in the plane of LISA, the variations in 
the sensitivity are maximum for the combinations Y^ and Y^ 2 \ From the Figures (Q) 
and (0), we note that the generators Y^ have zero sensitivity at few values of 0. This 
implies that no single data combination, even if it is a linear combination of these, can 
give optimal sensitivity. To obtain best results, one needs to maximise the SNR over 
the linear combinations amY^ 1 '. 

We analyse two strategies for optimising the sensitivity: 

(a) We consider the cyclic permutations of Y^ and compute the maximum sensitivity 
using these three data combinations. 

(b) We take the maximum of the linear combinations amY^ where am are complex 
numbers. At each value of 0, a different linear combination of Y^ is optimal. Thus 
the am for which the SNR is maximised are, in general, functions of 0. 

The results of these analyses are described in Figure @. For the strategy (a), the 
improvement in sensitivity averaged over is ~ 49 %; while for the strategy (b), the 
improvement in sensitivity averaged over is ~ 59 %. In addition Figure (S) shows the 
plots of am as function of 0, which maximise the SNR. 

The strategy (a) does not give the best SNR as it maximises the response over only 
a set of three data combinations. On the other hand, strategy (b) maximises the SNR 
over the module and is therefore superior. It is interesting to note that the maximised 
SNR is constant as a function of and is the same as the SNR of Y^ 1 ' maximised over 
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Figure 6. Plot of log S, of the and Y^ 2 \ are displayed in (a) and (b) respectively, 
as a function of 9 and (f> for / = 1 mHz and SNR=5 over a observation period of one 
year. 




Figure 7. Plot of sensitivity as a function angle <j) at 6 = ? and / = 1 mffz after 
averaging over the polarisations. 
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0. We remark that the strategy (b) gives only marginal improvement over strategy (a), 
and also strategy (a) is easier to implement requiring relatively fewer computations. 

As seen for a GW source with a frequency of 1 mHz, strategy (b) gives a constant 
value for the maximum sensitivity over (see Figure @). Thus the maximum sensitivity 
is independent of <fi for a given frequency. We further extend this analysis to other 
frequencies over the band-width of LISA, and obtain analogous results to the case of 
1 mHz. This maximum sensitivity however depends on the frequency. In Figure fllpp . 
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Figure 8. The coefficients a(j) maximising the SNR of the combination a^Y^ are 
plotted as functions of <j) for the frequency 1 mHz and 9 = ir/2. 
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Figure 9. Sensitivity curves for the toy model of LISA are given as functions of 
at / = \ raHz. The curves represented by broken lines correspond to the Y^ and 
its cyclic permutations. The solid line curves represent the sensitivities for the two 
strategies: (a) taking the maximum of Y^> and its cyclic permutations (thin solid 
line); (b) taking maximum over the linear combinations of Y^\ I = 1, 2, 3 (thick solid 
line). 
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Figure 10. Plot of sensitivity S as a function of frequency, for the toy model. 
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we present the maximum sensitivity curves as functions of frequency. For comparison, 
we perform analogous computations for the Michelson data combination X and the 
maximum of the cyclic permutations of In all these cases, the GW amplitude is 

first averaged over the polarisations and then for a fixed frequency, the sensitivity as a 
function of <fi is computed and then averaged over 0. Finally, this exercise is carried out 
for frequencies across the LISA band-width. 

The 1 mHz case is representative of the low frequency regime < 3 mHz where the 
improvement in sensitivity is about 49 % for strategy (a) and 59 % for strategy (b). 
These improvements scale up at higher frequencies at say, ~ 15 mHz, to 56 % and 67 
% , respectively, for the two strategies (a) and (b). 

5. conclusion 

In this paper, we have employed our previously set up formalism from paper I for two 
important applications: 

The first application is the extremisation of SNR. We have considered binaries as 
GW sources whose frequencies change at most adiabatically (monochromatic sources 
are included), that is, even if the frequency changes during the observation time, the 
change in SNRs under consideration is insignificant. The signal covariance matrix has 
been computed for which the signal is averaged over polarisations and directions. We 
have then computed the noise covariance matrix for LISA and obtained its eigenvectors 
and eigenvalues. This matrix has the same eigenvectors, resulting in calculational 
simplification. We have shown that the SNR for any data combination in the module, 



Max[a, Y 1 ] 
Michelson X 
Max[cyclic(Y 1 )] 
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lies between the upper and lower bounds which are determined by the eigenvectors 
of both matrices. We have further shown that the bounding SNR curves of the 
eigenvectors have multiple intersections within the band-width of LISA - 10~ 4 - 1 Hz. 
We have obtained the following results for the improvement in SNR: The improvement 
of SNR of the upper-bound over the Michelson combination goes upto 70 % , at 
high frequencies > 5 mHz, however, at low frequencies < 5 mHz, both have the 
same sensitivity. Since the eigenvectors are independent random variables, a 'network' 
SNR of independent detectors f| has been constructed from likelihood considerations 
which gives an improvement between a/2 and a/3 over the maximum of SNRs of the 
eigenvectors. The improvement of the network SNR over the Michelson combination is 
about 40 % at low frequencies < 3 mHz and rises above 100 % at high frequencies. 

The second application is a simple toy model of LISA rotating in its own plane. 
For this model we estimated the improvement of SNR by optimally switching from one 
data combination to another. We have shown that, if this strategy is used, it is possible 
to improve the SNR by about 55 % on an average over the band-width of LISA. We also 
show that if we instead maximise over the module, the SNR improves by about 60 % on 
an average over the LISA band-width. These improvements are obtained using just one 
data combination, namely, the eigenvector The SNR improvement in both cases 

is larger at higher frequencies > 10 mHz than at low frequencies. 
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